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A DIFFUSION  PHENOMENON  IN  ROTATING  AND  STRATIFIED  FLUIDS 

Sir  Geoffrey  Taylor,  in  one  of  his  first  papers  on  rotating  fluids 
(1921),  discovered  the  following  remarkable  phenomenon,  apparently  by 
accident:  Tf  a body  of  water  in  solid  rotation  is  disturbed  very  slightly 

and  a quantity  of  ink  of  the  same  density  is  dropped  into  it,  the  ink, 
instead  of  diffusing  in  a normal  manner,  gathers  into  vertical  walls. * 

The  walls  become  very  thin  and  very  elongated  as  time  goes  on,  until  they 
can  hardly  be  seen  when  viewed  with  a line  of  sight  perpendicular  to  the 
axis  of  rotation.  Viewed  from  above,  however,  they  have  the  appearance 
of  very  thin  ink  filaments  of  great  total  length. 

The  experiment  is  reproduced  in  fig.  1.  Just  before  the  ink  was 
inserted,  a weak  relative  motion  was  generated  by  a gentle  stirring. 

The  Rossby  Number,  using  the  radius  of  the  pan,  is  about  0.03.  The 
angular  velocity  is  10  rpm.  For  the  purposes  of  comparison  a similar 
experiment  was  performed  in  a non-rotating  body  of  water.  The  ordinary 
diffusion  process  is  shown  in  fig.  2. 

Taylor's  explanation  of  this  phenomenon  was  incomplete,  and  the 
purpose  of  this  note  is  to  advance  a qualitative  explanation  which  seems 
to  serve  the  purpose.  Taylor  noted  that  the  relative  velocities  were  so 
small,  compared  to  the  basic  rotation,  that  steady  motions  had  to  be  two- 
dimensional  , i.e„,  no  vertical  motion  and  no  vertical  shear.  This  follows 

^It  is  not  necessary  to  disturb  the  water  before  putting  in  the  ink. 

The  walls  will  form  without  the  added  disturbance,  but  they  remain  so 
close  together  that  they  are  more  difficult  to  see  in  a photograph. 


(e)  1 min.  12  sec.  (f)  5 min. 

Figure  1.  Taylor's  ink  walls. 


(e)  48  sec. 


(£)  5 min. 

Figure  2.  Diffusion  in  a non-rotating  liquid. 


- 4 - 


from  the  steady  state  vorticity  equations  of  a homogeneous  frictionless 
rotating  fluid: 

ufx+ vCy  + w£z  s £ux  + ijuy+  (f  +Uuz 
Ul?j  + Vl?y+  Wl?2  s fvx+,vy+(f  + Uv2 
utx+ v£y+ w£z  = f wx+7?wy+(f+^)w2 

where  the  relative  vorticity  is  ) and  the  axis  of  rotation  is 

parallel  to  the  z-axis.  As  the  relative  velocities  become  smaller  and 

smaller  all  terms  involve  products  of  small  quantities  except  fuz,  fv_, 

fw_.  If  we  assume  that  derivatives  of  velocity  do  not  remain  finite  as 

the  velocities  become  infinitesimal,  u , v , and  w must  become  second- 

' z*  z9  z 

order,  and  ultimately  negligible  compared  to  derivatives  of  velocities 

I 

in  the  horizontal  direction.  Furthermore  in  our  experiment  w is  zero 
at  the  bottom  of  the  vessel.  Since  its  vertical  variation  is  negli- 
gible, w is  zero  everywhere. 

The  above  argument  shows  that  once  a vertical  wall  of  ink  is 
formed  it  will  remain  as  such,  undeformed  by  subsequent  motions.  For 
the  time  intervals  of  the  experiment,  molecular  diffusion  is  nearly 
negligible. 

Taylor’s  explanation  is  quite  correct,  but  does  not  explain  why 
the  ink  gathers  originally  into  walls,  nor  why  the  walls  become  so 
thin  and  elongated  vertically  and  horizontally.  A rotating  fluid 
possesses  vorticity  and  this  is  represented  by  vortex  lines  parallel 
to  the  axis  of  rotation  if  there  is  no  relative  motion.  The  number  of 
lines  per  unit  horizontal  area  is  a measure  of  the  intensity  of  rota- 
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tion.  Neglecting  friction,  vortex  lines  are  material  lines;  in  all 
subsequent  motions  they  are  composed  of  the  same  fluid  particles.  If 
we  suddenly  plunge  a foreign  body  of  fluid  into  a system  in  solid 
cyclonic  rotation,  the  vortex  lines  will  be  forced  apart,  and  as  a 
result  the  relative  circulation  on  horizontal  rings  of  fluid  surround- 
ing the  foreign  fluid  will  become  negative.  Initially,  the  resulting 
anticyclonic  vortex  will  not  have  a balance  of  pressure  gradient  and 
Coriolis  force,  and  the  rings  will  contract.  This  will  force  the 
injected  fluid  to  spread  vertically  until  it  reaches  the  bottom  of 
the  vessel  and  the  free  surface.  The  relative  motion  will  now  be 
two-dimensional  and  the  newly-formed  ink  wall  will  be  stretched  out 
horizontally  by  local  shearing  motion. 

Such  a remarkable  type  of  diffusion  process  might  seem  at  first 
glance  to  be  a most  unusual  occurrence,  or,  indeed,  peculiar  to  rota- 
ting fluids.  This  is  not  the  case,  however;  in  fact,  the  above  explan- 
ation occurred  to  me  while  reflecting  on  the  tendency  for  cigarette 
smoke  in  a quiet  room  to  concentrate  into  sheets  of  smoke  parallel  to 
the  floor.  A similar  phenomenon  exists  in  the  atmosphere  when  smoke 
is  discharged  into  the  air  during  stable  conditions.  Using  the  above 
approach,  the  injection  of  the  smoke  causes  a local  separation  along 
the  vertical  of  the  surfaces  of  constant  density.  The  action  of  grav- 
ity, in  attempting  to  restore  the  horizontal ity  of  these  surfaces,  will 
compress  the  smoke  vertically  and  spread  it  horizontally  into  sheets. 
The  analogy  to  the  rotating  phenomenon  is  very  close. 
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A more  intuitive  explanation  of  the  formation  of  sheets  of  smoke 
in  a stable  atmosphere  is  that  the  smoke  rises  under  the  action  of 
gravity  until  it  reaches  an  inversion  surface.  This  prevents  any 
further  vertical  motion  and  the  smoke  spreads  horizontally.  In  order 
to  check  the  analogy  to  a rotating  system,  and  to  dismiss  the  neces- 
sity for  a near-discontinuity  in  density  gradient,  implied  by  the  last 
argument,  an  experiment  was  conducted  with  ink  injected  into  a stable 
water  and  salt  solution.  By  careful  preparation  a vessel  was  filled 
with  a mixture  having  a smooth  density-height  curve  (fig.  3).  A 
quantity  of  ink,  with  a density  equal  to  the  mean  density  of  the  fluid 
in  the  vessel,  was  then  injected.  The  subsequent  developments  are 
shown  in  fig.  4.  The  ink  forms  a fairly  thin  sheet  within  approxi- 
mately 5 min.  After  a half-hour  or  so,  the  gravitational  diffusion 
effect,  tending  to  spread  the  ink  horizontally  and  contract  it  verti- 
cally, is  reduced  to  the  level  of  importance  of  molecular  diffusion. 
Before  this  occurs,  the  diffusion  process,  if  it  were  Fickian, 

(Sutton,  1953),  would  involve  large  horizontal  diffusion  coefficients 
and  a negative  vertical  coefficient. 

It  is  of  some  interest  to  investigate  the  time  interval  required 
for  the  formation  of  the  ink  walls  in  a rotating  fluid.  This  may  be 
done  crudely  by  a dimensional  argument.  Assuming  that  a spherical 
volume  V of  foreign  material  is  injected  into  an  unlimited  rotating 
fluid  of  the  same  density,  and  denoting  the  decreasing  horizontal 
dimension  by  d,  we  have  simply 

d=  v3x(at). 


DEPTH  (inches) 
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Figure  3.  Density-height  curve  for  experiments  of  fig.  4. 


Figure  4.  Formation  of  ink  sheets  in  a stratified  fluid 


Si  i.  m 
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This  assumes,,  of  course,,  that  the  volume  retains  its  initial  axial 

symmetry  so  that  it  spreads  vertically  and.  not  horizontally.  The 

-■) 

experiments^"  indicate  that  d changes  by  a factor  of  1/10  in  perhaps 

-I  "5 

30  sec  at  an  angular  velocity,,  XI  - • i sec~  . Thus  a change  to  d/V 

of  1/10  is  associated  with  a value  of  30  for  XI  t,  In  the  atmos- 

-5  - ] 

phere,  where  XI  is  about  7x10  sec  , such  a change  would  require 
t = 30/7 xlO"5 sec  or  about  5 days.  For  any  conceivable  atmospheric 
phenomenon  of  this  kind,  the  effect  of  eddy  diffusivity  would  be 
several  orders  of  magnitude  greater. 

The  analogous  effect  in  a stratified  fluid  would  seem  to  be  of 
more  practical  importance.  A dimensional  analysis  of  this  case  leads 
to 


d = 


I 

V*a 


.2 

gt  a \ 

Po  ' 


» 


where  d is  the  diminishing  vertical  dimension,  pa  is  the  uniform 
density  of  the  foreign  matter,  a is  the  density  (potential  temper- 
ature) gradient  of  the  environment,  and  g is  gravity.  The  first  term 
in  the  argument  of  <f>  is  the  non-dimensional  density  difference  be- 
tween the  ink  mass  and  its  environment.  If  we  regard  this  quanti  ty 
as  important  only  insofar  as  it  is  associated  with  the  buoyancy 
effect  (and  therefore  only  when  associated  with  g),  we  may  neglect 
this  number.  In  the  experiment  a/pa  was  about  3xl0'*^cm“^,  and  the 

“The  experiment  is  limited  vertically  by  the  bottom  of  the  vessel 
and  the  free  surface.  A full  analysis  should  include  this  dimen- 
sion and  gravity. 
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time  t for  a reduction  of  d to  1/10  of  its  original  value  was  a'uout 
5 min.  If  we  choose  a mass  in  a stable  environment  with  a potential 
temperature  rise  of  5-1CC  in  a hundred  meters,  the  requirement  that 
the  second  argument  of  $ be  equal  in  the  laboratory  and  the  atmos- 
phere yields  a value  of  1-3  hours  for  the  mass  in  the  atmosphere  to 
form  a sheet.  This  assumes,  of  course,  a complete  absence  of  eddy 
diffusivity. 

This  time  period  seems  toe  long  if  the  above  process  is  respon- 
sible for  the  frequent  tendency  for  smoke  to  concentrate  into  sheets 
in  a stable  atmosphere.  R.  ports  of  this  phenomenon  (Barad,  1951) 
indicate  that  it  occurs  in  a matter  of  minutes.  It  is  possible 
that  the  local  stability  at  the  height  of  the  stack  may  be  consider- 
ably greater  than  that  assumed  in  the  above  calculation.  It  is  more 
likely,  however,  that  the  above  analysis  is  too  crude.  In  any  event, 
the  vertical  eddy  diffusivity  drops  to  exceedingly  small  values  if 
the  stability  is  high.  It  is,  therefore,  reasonable  to  assume  that 
the  gravitation  effect,  tending  to  condense  vertically  and  spread 
horizontally,  is  frequently  of  the  order  of  magnitude  or  greater 
than  the  vertical  eddy  diffusivity.  This  effect,  may,  moreover,  help 
to  explain  the  findings  by  Parr  (1936)  that  stability  increases  hori- 
zontal diffusivity. 
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APPENDIX 

It  seems  desirable  to  present  a theoretical  development  which  helps 
to  explain  the  analogy  between  rotating  and  stratified  fluids  with  re- 
spect to  the  diffusion  phenomenon  described  above.  We  divide  the  inves- 
tigation into  two  parts: 

A.  Stratified  fluid. 

According  to  the  description  of  the  experiments  the  introduction  of 
ink  serves  initially  to  distort  the  density  distribution  in  one  case, 
and  the  pattern  of  vortex  lines  in  the  other.  In  this  portion  of  the 
investigation  we  wi 1 1 assume  a disturbance,  p',  of  the  density  field. 

The  latter  is  given  by  pQ  e”^z,  when  undisturbed.  The  coordinates  are 
z,  upward,  x,  horizontal.  For  simplicity  we  assume  that  the  disturbance 
is  independent  of  the  y-direction.  This  does  not  conform  to  the  situation 
in  the  experiment  or  in  the  atmosphere  but  will  serve  sufficiently  well 
for  the  present  purpose.  The  exponential  density  distribution  was  chosen 
for  simplicity.  If  we  take  the  origin  of  the  z-axis  near  the  disturbance, 
0z  will  be  very  small  in  the  vicinity  of  the  disturbance  and  the  density- 
gradient  will  be  nearly  linear  in  this  region.  In  fact,  we  will  assign 
a length  L for  the  characteristic  vertical  dimension  of  the  disturbance 
and  assume  /9L«1.  Then,  if  zq  is  the  undisturbed  height  of  a given 
density  surface,  p'  = po(e”^zo  - e“^z)  — pQ/9(z  - zQ)  = o(pQ/9L).  From 
this  we  infer  that  p*/po«l,  an  entirely  reasonable  result  of  the  assump- 
tion, /JL«1. 


The  vorticity  equation  is 


iL  + Xp  jLb. _ _L p 

dt  e * dt  p z dt 


C = »x-u 


z 


t 


(1) 

(2) 


where  partial  derivatives  are  denoted  by  subscripts  as  before.  Since  the 
pressure  distribution  is  very  closely  hydrostatic  we  may  neglect  the  sec- 
ond term,  compared  to  the  last  term  in  equation  (1).  We  assume  now  that 
the  motion  is  very  slow  and  neglect  the  non-linear  terms  in  the  velocity 
and  density  perturbation.  Equation  (1)  becomes 


a t p,  p* 

dt  P \ P 


g = 0 . 


(3) 


By  the  same  assumptions  the  requirement  of  incompressibility  leads  to 

P = w />  /3  . M) 

t o 

Introducing  the  stream  function  by  the  expressions 

u = , w=*  (5) 

4 ' A t 


the  vorticity  equation  becomes 


vV?fztl+  = *zx+  *zz  ■ 


(6) 


The  ratio  of  the  second  term  to  the  first  term  is  of  the  order  of  /3L,  so 
that  we  have,  approximately. 
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Since  initial  conditions  are  given  in  terms  of  the  disturbed  positions 
of  the  density  surfaces,  we  may  introduce  the  quantity  j^=  z - zQ.  We 
obtain  finally 


B.  Rotating  fluid. 

Assuming  axial  symmetry  of  the  disturbance  about  the  x-axis  (axis 
of  rotation),  the  equations  of  motion  are 


d v 
dt 


*--2Slw 


dx 

fir  * 


(9) 


du_  . fix 

dt  6>x  » 


(10) 


i.(flr2+wr)  _ n , 


(11) 


where  X=  p/p  plus  centrifugal  and  external  forces,  v is  the  radial  ve- 
locity, u is  the  velocity  along  the  x-axis,  w is  the  velocity  tangent  to 
circles  about  the  x-axis,  and  r is  radial  distance. 

The  vorticity  equation  is 


r 


2w  Wy 

r 


-2ilwx  = 0 # 


(12) 


f 


vx” 


(13) 
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Neglecting  second-order  terms  in  the  velocities,  this  becomes 

£t-2ftwx=C,  (14) 

and,  to  the  same  order  of  approximation,  (11)  reduces  to 

wt+2ftv  = 0.  (15) 

Introducing  the  stream  function  by  the  expressions 

u = - *r,  v=  (16) 

and  eliminating  w,  the  vorticity  equation  is 

v2+=+kx+V'7L  • <17) 

We  now  introduce  8‘  = r^  - rQ~,  where  ro  is  the  Lagrangian  distance  of 
a vortex  line  from  the  x-axis.  This  expression  occurs  in  the  integral 
of  (11): 

ilr2+  wr  = . (in) 

In  terms  of  8*,  (17)  becomes 

v2*ttt+4“2sm=o  <19) 

It  is  remarkable  that  equations  (8)  and  (19),  governing  the  develop- 
ment of  an  initial  disturbance  in  a stratified  and  rotating  fluid,  respec- 
tively, are  almost  identical  except  for  the  fact  that  the  Laplacean  opera- 
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tor  is  in  cartesian  coordinates  in  one  case  and  in  cylindrical  coordi- 

O 

nates  in  the  other.  The  parameters  g£  and  4X1“  show  that  tiie  time  scales 
of  the  two  phenomena  are  in  the  ratio  ^ 

Ts/Tr=  Za/igP)2  (20) 

so  that  geometrically  similar  initial  disturbances  should  develop  in  time 
periods  having  this  ratio,  provided  the  fluids  are  effectively  infinite  or 
the  boundary  conditions  are  geometrically  similar.  For  the  experimental 
conditions  of  this  paper  this  ratio  is  about  one,  whereas  the  experimental 
results  indicate  roughly  that  the  ratio  is  about  10  (300  sec/30  sec).  Pos- 
sibly this  results  from  the  fact  that  the  ink  was  injected  along  a line 
perpendicular  to  the  density  surfaces  in  one  case  and  parallel  to  the  vor- 
tex lines  in  the  other. * It  is  physically  clear  that  this  would  make  a 
considerable  difference  in  the  experimental  results.  The  theoretical  in- 
dication that  the  time  period  in  the  stratified  case  may  be  as  little  as 
1/10  of  the  1-3  hr  period  computed  from  dimensional  considerations  of  the 
experimental  results,  is  interesting  with  regard  to  smoke  diffusion  in  the 
atmosphere.  If  smoke  issues  from  a stack  at  a temperature  close  to  that 
of  the  air,  the  injection  would  be  more  or  less  along  a line  parallel  to 
the  density  surface  and  the  vertical  flattening  may  occur  in  a matter  of 
minutes  instead  of  the  1-3  hr  period  suggested  by  the  experimental  evidence, 

I’Ve  assume  here  that  the  different  meanings  of  Sand  8'  (z  - zQ,  r~  - r()~), 
and  the  two-dimensionality  of  the  stratified  portion  of  the  theory  are  in- 
sufficient to  explain  the  order-of-magnitude  difference  of  the  time-scale 
ratios  in  experiment  and  theory,  A more  involved  theoretical  investigation 
would  be  needed  to  decide  this. 
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